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Abstract 

A nonequilibrium statistical operator method is developed for ensembles of particles obeying 
non-Hamiltonian equations of motion in classical phase space. The main consequences of non-zero 
compressibility of phase space are examined in terms of time-dependent dynamic quantities. The 
generalized transport equations involve the phase-space compressibility in a non-trivial way. Our 
results are useful in molecular dynamics simulation studies as well as nonequilibrium or quasiclas- 
sical approximations of quantum-classical dynamics. 
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Theoretical interest in the non-Hamiltonian motion appears in the early works 
where it was realized that the constant-temperature simulations require the introduction of 
additional non-Hamiltonian acceleration into dynamical equations. This has opened new re- 
search directions associated with extended systems. For the past few years the physics of non- 
Hamiltonian systems was partially developed for molecular dynamics simulation purposes 
4], 0, 0, Q, . Particular applications of non-Hamiltonian dynamical systems are overviewed 
in a recent paper 9( , where the authors exploit their formulation of non-Hamiltonian statis- 
tical mechanics 101. In this regard it is important to mention also the interesting studies of 
mlX e d q — classical SyStemS fiflQ, where tbe aca-Haaultcaiaa behavicr 

plays an essential role. 

One of the most popular non-Hamiltonian ensembles is the so-called Nose- Hoover chain 
In case with two thermostat variables, the corresponding Hamiltonian is 
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where rji and p m are the conjugate thermostat variables, g is the number of degrees of 
freedom, T is the thermostat temperature. It can be shown that the Hamiltonian ((TJ) is 
conserved by the following (non-Hamiltonian) equations of motion: 
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In real MD simulations the particular choice of non-Hamiltonian ensemble can vary de- 
pending on the physical problem under study. But all of them can be described in general 
formalism we present here. Recently using a generalization of the symplectic form of the 
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Hamilton equations of motion the authors of Ref. Q showed that there is a unique gen- 
eral structure that underlies the majority of the equations of motion for extended systems. 



Furthermore, Sergi 17j, [l8| proposed a consistent theory of statistical mechanics for non- 
iamiltonian systems under equilibrium conditions. The tempting feature of the theory 
lfl [nl . Q| is that it is based on the introduction of a suitable Poisson bracket, and has 



the general structure of classical formalism 
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with a natural appearance of the nonzero 



phase-space compressibility. The main feature of the systems of interest is energy conserva- 
tion, which is a reason of non-Hamiltonian behavior. It is interesting therefore to consider 
the situation, when the energy is still conserved, but the system itself is out of equilibrium, 
particularly when the dynamical quantities depend on time explicitly. 

The starting point in describing the non-Hamiltonian statistical ensembles is the intro- 
duction of the proper non-Hamiltonian Poisson bracket: 

ik 

where x = (q,p) = (qi, <2 , 2---<?v; Pi,P2---Pn) is the 2iV-dimensional vector in phase space, A(t) 
and C(t) are arbitrary dynamical quantities (we allow them to depend on time explicitly), 
Bik = is an antisymmetric matrix, whose elements are functions of x and, possibly, 

time. Following Ref. [l7| we write the equations of motion in standard form ij = {xj,Ti.}, 
where 7i is the Hamiltonian, which is a constant of motion by construction due to of anti- 
symmetry of the matrix Bik- 

The main feature of the non-Hamiltonian system is the failure of the Jacobi relation, 
which implies time-irreversibility of non-Hamiltonian mechanics. This is formally a conse- 
quence of our requirement for eq. (JZJ) and the standard equation of motion to be satisfied 
simultaneously. By defining the Liouville operator 

f A „ dAOH fAnn 



we obtain the generalized Liouville equation (see Ref. Il6l ) in the form 



jLp(q,P;t) = - (iL + p(q,p;t), (8) 
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^ dxi dxj 



is the phase space compressibility. For the purpose of generality we also allow it to depend on 
time explicitly (through Bij). For particular determination of the matrix Bij (and therefore 
x) it is useful to write the equation of motion in the form 
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For the particular ensemble, described by Eqs. ((TJ - (JHJ), we obtain 
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so that the compressibility can be calculated as 
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In general for a system of N identical particles (more exactly for N degrees of freedom) 
we define the statistical time-dependent average as 



( A ) t = ^2 / A N(q,P;t)p N (q,p;t)dT N = Ti {A(q,p;t)p(q,p;t)}, 



(9) 



where Apf(q,p;t) denotes the dynamic variable A for a system of N particles, i.e., 
A(qi, ...,q N ;p 1 , ...,p N ;t), dT N is the element of phase volume. 

The time evolution of an arbitrary dynamical variable obeys the usual equation 



dA dA r . nn dA . f . 
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One interesting consequence of the time-irreversibility is that the average derivative of a 
dynamical variable is not necessarily equal to the derivative of its average. In fact, taking 
into account eq. (JBJ) we obtain the relation 

d(A) t 
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Equations and (|12|) are particularly important in identifying the real measured mean 
quantity. Another important peculiarity of non-Hamiltonian systems follows if one is in- 
terested in the time evolution operator for the initial statistical distribution p(q, p; t) = 
U^(t, t f )p(q,p; 0). Following the standard procedure (^(J we obtain the time-ordered evo- 
lution operator U± (t, t') in the following form 

U^\t,t') = exp ± I- J (iL(r) + x) drj (13) 

where exp ± denote the correspondingly time-ordered exponents. We emphasize that the 
evolution operator in Eq. ()13j) does not describe the evolution of dynamic variables. Eq. 
(PHI) gives A(q,p;t) = Ul(t,t')A(q,p;t'), where 

Ul{t,t') = U^ ] \t,t') = exp ± li J L(r)drJ . (14) 

Let us suppose that the given nonequilibrium non-Hamiltonian ensemble can be char- 
acterized by the average values (Cn)* °f a finite set of dynamical variables {( n {t)}- Then 
the time evolution of the nonequilibrium macroscopic state is described by the generalized 
transport equations (see Eqs. (fTT|) and (fT2l ) 

d(CnY 



dt 



Tr { ((„ - xC„) p(t)} = Tr {({(„, H} - p(t)} , (15) 



where the trace operation is defined in Eq. ©. Let p(t) be a trial statistical operator, such 
that (Cq) 1 = Tr (Cnp(^)), and Tr (p(t)) = 1. We choose the entropy functional in the standard 
form 20j . 

S[p(t)] = -Tr (p(t)lnp(t)) -Y,<Pn(t)Tr(( n p(t)) ~ W^ffl, (16) 

n 

where <p n {t) and X(t) are the corresponding Lagrange multipliers. Varying (J 16)) with respect 
to the trial distribution p(t) we obtain the relevant distribution, 

Pr=exp j-$(t)-^ n (t)C„j, (17) 



where is t 
Zubarev 20, bj 



re Massieu-Plank function. Here we use the notion "relevant" following 
for the distribution that maximizes the entropy functional ()16|) at a 



particular time. The relevant distribution is not yet the required nonequilibrium distribution 



because, in general, it does not satisfy the Liouville equation. Nevertheless, as shown below 
the relevant statistical distribution will serve as an auxiliary distribution to determine special 
solutions of the Liouville equation that describe irreversible macroscopic processes. 
The Lagrange multipliers (p n are determined from the self-consistency conditions, 

(CnY = (CnYr = K (CnPr(t)) , (18) 

so that S(t) = $(£) + J2 n <Pn(i) (CnY- The meaning of the noneq uilibrium entropy is the 
same as the corresponding entropy in the Hamiltonian systems (20, |2l| . In contrast with 
the Gibbs entropy, S(t) can change in time. Let us demonstrate the specific feature of this 
entropy in the non-Hamiltonian system by calculating the derivative 

ds(t) ^ 5s(t) d(( n y ^ d^y ^ r/z t \ M \ (m , 

n N ^ n/ n n 

It seems confusing that dS(t)/dt ^ even for constant ( n due to the compressibility of the 
phase space. This of course is misunderstanding, because the entropy and its evolution are 
essentially determined by the mean values of ( n and by the derivatives of their mean values. 
This is particularly important in applications of the theory since the interpretations of the 
desired quantities depend on the definition of the averaging procedure. 

To find the formal solution of the generalized Liouville equation (jHJ), we assume that at 
some initial instant of time the distribution coincides with the relevant one (|17|1. satisfying 
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the self-consistency relations (jl8j) . Following the method proposed by Zubarev 
retarded solutions of the Liouville equation in our case is found by solving the following 
perturbed equation 

^ + (it + x) p{t) = -e (p(t) - p r (t)) (20) 

with e — > in the final expressions. The formal solution of the initial Liouville equation (jHJ) 
reads 

i 

p(t)=\im Q e J e-^^e-^-^ 1 p r {t')dt' . (21) 

— oo 

Integrating eq. (j7T|) by parts we arrive at 

t 

p{t) = p r (t) - lfme J e -(e+x)(t-t') e -i(t-t')L + X ^J pr {t')dt'. (22) 
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In principle equation (}22j) answers the question how to construct nonequilibrium distribution, 
which is a functional of dynamic variables and the phase space compressibility, if the system 
develops from a relevant state. 

Let us further write the statistical operator as a sum of relevant distribution and local- 
time deviation pit) = p r {t) + Ap(t), so that equation (|2U|) becomes 

^ + it + x + ej Ap(t) = - ^ + it + xj p r (t). (23) 

Then, keeping in mind that p r (t) depends on time only through the relevant variables ( n , 
we can calculate its time derivative 

dp r (t) \- 6p r (t) d(( n y ^ Sprjt) {{■ \ 



dt ^6(( n y dt ^5(CnY 

-E^fcT}-?^^^^)^}- (24) 

We define the Kawasaki-Gunton operator [23] in the standard manner 

V KG A = Pr (t)TrA + 6 Sf t (Tr(AC„) - (Q^A} . (25) 

Using the relations Tr{iLp(t)} = 0, Tr{p(t)} = — (x)', one can prove that the projection 
operator (21) conserves its features in the non-Hamiltonian case, namely VKG{t)p{t) = p r (t), 
T :> KG(t)p(t) = p r {t). Equation (j2U) can now be rewritten in the projected form 

^ = -V KG (t)(iL + K) P (t). (26) 

Inserting equation (|2l)j) into Eq. (|23|) and integrating over time, we obtain a solution for 
Ap(t) as follows 

t 

Ap(t) = - J e-^U^\t,t')Q KG {t') (it + x) p r (t')dt', (27) 

— oo 

where Qkg = 1 _ Vkg-> and (compare with Eq. (fTTfy ) 

U& (t, t') = exp | - J Q KG {t) (it + x) dr I . (28) 



Equation (|27|) allows us to express the right-hand side of Eq. (|15|) in terms of relevant 
distribution, i.e. in terms of variables (Cn)*- So we get the closed set of generalized transport 
equations for observables. 



To obtain more transparent equations we define the entropy operator 



so that p r {t) = exp{— Introducing the Mori projection operator 24 1 



V M (t)A = (A)l + J2 |^ (Cn - (&>') , (29) 
which satisfy the identities V^it) = "Pm(^), Vm^Cu — Cn, one can prove the relation 

Q KG (t) {(iL + x)p r } = Q M {t) jx - iLs} pr, (30) 

where Qm(^) = 1 — 'Phdif) is complementary to the Mori projector (j29j). Now the generalized 
transport equation ([To]) reads 

di {Cn)t = di {Cn)tr + ^ { Qu{t) " xCn ) Ap ) ' (31) 
where Ap must be taken from equation (J27|) . It is useful to define new dynamical variables 
(generalized fluxes) 

T n {t) = Quit) (Cn - <„) , 

which are orthogonal to the space of relevant variables in a sense that VuT-rSfy — 0. The 
variables T n (t) determine the effects of macroscopic degrees of freedom of the evolution 
on the nonequilibrium state of the non-Hamiltonian system described by a set of relevant 
variables and a given phase space compressibility x. 

Incorporating equations ([27)1 . (j2UJ), and (J3T|) we obtain the final result for the generalized 
transport equations 

t t 
■(Cn) i = ^(Qt + E / e^-^^M'K^' + / e-^V^M')^', (32) 



dt 



where 

V mn (t,t') = Tr {l m (t)U^(t,t')l m (t')p r (t')} (33) 
are the normal generalized kinetic coefficients, and 

FmnM = Tr ^I m (t)U^\t,t')Q M (t')H (y^ifnCn ~ Pr(0 j ( 34 ) 



are the effective forces determining the anomalous contribution into the system evolution 
associated with nonlinear phase space compressibility. 

The last term in Eq. (J32j) is associated with the initial irreversibility of the process due 
to non-Hamiltonian behavior, while the second term is due to the nonequilibrium process 
itself. It does not vanish even when x = 0, and coincides with the well-known result for 
Hamiltonian ensembles. These terms demonstrate the general structure of the solution and 
clearly show the contributions from two physically different processes. In practice equations 
(I33|) and ()34|) can be used with standard simplifications for particular systems. 

In conclusion, we have proposed a possible realization of the nonequilibrium statistical 
mechanics for ensembles of particles obeying non-Hamiltonian dynamics. Starting with the 
non-Hamiltonian Poisson bracket (J7|) and the equations of motion in standard form, we 
derived the basic relations for time-dependent dynamical variables with conserved Hamilto- 
nian. Introducing the reduced description method in analogy with Hamiltonian dynamics 
we found the retarded solution (}22|) of the generalized Liouville equation (JEJ). We emphasize 
however, that Eq. (f2"2"J) is the specific solution, not general, and it describes the evolution 
from a specific auxiliary relevant distribution (this is reflected in the introduction of pa- 
rameter e). Using the projection operator technique we obtained the generalized transport 
equations (}3"2j) . which contain two physically different contributions due to the statistical ir- 
reversibility ()33|) and the time-irreversibility (}3"4")) associated with non-Hamiltonian features 
of the system. 
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